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We investigate the previously proposed cyclic regime of the Kosterlitz-Thouless renormalization 
group (RG) flows. The period of one cycle is computed in terms of the RG invariant. Using 
bosonization, we show that the theory has U q (sl(2)) quantum affine symmetry, with q real. Based 
on this symmetry, we study two possible S-matrices for the theory, differing only by overall scalar 
factors. We argue that one S-matrix corresponds to a continuum limit of the XXZ spin chain in 
the anti- ferromagnetic domain A < — 1. The latter S-matrix has a periodicity in energy consistent 
with the cyclicity of the RG. We conjecture that this S-matrix describes the cyclic regime of the 
Kosterlitz-Thouless flows. The other S-matrix we investigate is an analytic continuation of the 
usual sine-Gordon one. It has an infinite number of resonances with masses that have a Russian 

' doll scaling behavior that is also consistent with the period of the RG cycles computed from the 

beta-function. Closure of the bootstrap for this S-matrix leads to an infinite number of particles of 

. higher spin with a mass formula suggestive of a string theory. 
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^ . I. INTRODUCTION 

CO ' 

^ ■ The Renormalization Group (RG) continues to be one of the most important tools for studying the qualitative and 
quantitative properties of quantum field theories and many-body problems in condensed matter physics. A widely 
important class of theories are those with RG fixed points in the ultraviolet (UV) or infrared (IR) , and our intuitive 
understanding of the generic behavior of quantum field theory is largely based on theories with these properties. The 
notion that massive states decouple in the flow toward the IR is an example of such an intuition. However fixed point 
behavior is not the only possibility, and physically sound examples with other kinds of behavior are important to 
explore. This paper is concerned with cyclic RG flows, whose possibility was considered as early as 1971 by Wilsonp]. 
~- -. However at the time no interesting models were known that exhibited this behavior. 

Recently a cyclic RG behavior has been found in a number of models, wherein the couplings return to their initial 
^ i values after a finite RG time A: 

g(e x L)=g(L), (1) 

where L is the RG length scale 0, 0i QL ■ The models in 0, Q are problems in zero-dimensional quantum mechanics. 
The model in is a natural extension of the BCS model of superconductivity and thus a many-body problem. 

The cyclic RG property eq. Q has some important implications for the spectrum of the hamiltonian. Namely, if 
{E n , g, L} is the spectrum of eigenvalues of the hamiltonian for a system of size L, then 

{E n , g, e x L} = {E n , g,L}, (2) 

i.e. the energy spectrum at fixed g should reveal a periodicity as a function of L. A nice feature of the model in [B| 
is that in addition to the beta-function with the cyclic property one could obtain analytic results for the spectrum 
using the standard BCS mean field treatment. In this way one could study the interplay between the cycles of the 
RG flow and the spectrum. The manner in which the spectrum reproduces itself after one RG cycle is dependent on 
the existence of an infinite number of eigenstates with the 'Russian doll' scaling behavior in an appropriate limit 

E n +i ~ e x E n . (3) 

In each cycle these eigenstates reshuffle themselves such that the (n + l)'th state plays the same role as the ro'th state 
of the previous cycle. 
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FIG. 1: Kosterlitz-Thouless flows. 



Another possible signature of a cyclic RG is a cyclicity of the S-matrix itself. Applying RG equations to the S-matrix 
S(E), where E is the energy, one expects: 

S(e x E) = S(E) (4) 

The above equation was anticipated in p). 

The models considered in [j] are relativistic models of quantum field theory in 2d. What is surprising is that 
this theory is in fact a well-known theory that arises in a multitude of physical problems: an anisotropic left-right 
current-current interaction that gives rise to the Kosterlitz-Thouless (KT) flows. The standard picture of these flows, 
based on the one-loop approximation, is shown in figure 1. The g\\ axis is a line of IR or UV fixed points. The flows 
in the region \g±\ > \g^\ on the other hand are already peculiar in that they both originate and terminate at no 
clearly identifiable fixed point. In £| it was proposed that this region has a cyclic RG flow based on the all-orders 
beta- function conjectured in 6J. However since the flows extended beyond the perturbative domain j_| < 1, it 
remained unclear whether other non-perturbative effects would spoil the cyclicity. Furthermore, at the time little was 
known about the spectrum. In any case, if this proposal is correct it implies that an important feature of the KT 
flows has been overlooked. 

The main purpose of this article is to provide further evidence for the cyclicity of the KT flows by proposing an 
S-matrix. Let us summarize our main results. In theories with UV or IR fixed points the universal properties encoded 
in the beta-functions are anomalous dimensions of fields which are related to the slope of the beta-function at the 
fixed point. On the other hand, for a cyclic flow, the only universal, coordinate independent property is the period 
of one cycle. We compute the period A in section II in terms of the RG invariant Q = ~h 2 /16 (at small coupling 
Q w <7jj — <?^), and find the result A = 2w/h. 

In section III we bosonize the currents and obtain a sine-Gordon theory with a special form of the sine-Gordon 
coupling with both real and imaginary parts. (See eq. Other regimes of the KT flows are known to have an 

exact S-matrix description corresponding to the sinh-Gordon and massive/massless sine-Gordon theories. In section 

IV we argue that the theory possesses U q (sl(2)) symmetry with q real. (The usual sine-Gordon model corresponds to 
q a phase.) This symmetry fixes the S-matrix up to an overall scalar factor. The only free parameter of this S-matrix 
is h which governs the period of the RG cycles, and we find q = —e~ 7rh / 2 . 

For the simplest choice of the overall scalar factor, our S-matrix can be understood as a certain relativistic continuum 
limit of the XXZ spin chain with A < — 1. In terms of A, the period of a cycle is 

A= . (5) 

coBh _1 (-A) 

In this case, the signature of the cylic RG is the periodicity of the S-matrix eq. 0, where the period agrees precisely 
with the beta function computation in section II. 

The other possible overall scalar factor we investigate is an analytic continuation of the usual sine-Gordon one. The 
resulting S-matrix is not real analytic, and for this reason we believe the S-matrix related to the XXZ spin chain is the 
correct one physically. Nevertheless this second S-matrix has some interesting properties we explore in the last section 
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of the paper. The soliton S-matrix possesses poles which imply the existence of an infinite number of resonances with 
mass 

T17T 

m n = 2M S cosh — , n = 1, 2, 3, oo, (6) 

where M s is the mass of the soliton. We show that these resonances are actually a spin 1 triplet and a singlet. The 
above spectrum, eq. JHJ, which is rather novel in the subject of integrablc quantum field theory, is the Russian doll 
spectrum that was anticipated based on the analogy with the results in Q. Namely, 

m n ~ M s e n7r/h m n+2 ~e x m n , n > -, (7) 

to be compared with eq. J3J). This property of the spectrum is a strong indication of the cyclicity of the RG. 
Furthermore, we find that closing the S-matrix bootstrap leads to resonances of higher spin with a mass formula that 
suggests a string theory description. Theories with an infinite number of resonance poles were also found in 0, @ ; 
however these theories appear to be rather different, the S-matrices being built out of elliptic functions. 



II. CYCLIC RG FLOWS IN ANISOTROPIC SU{2) CURRENT PERTURBATIONS 

We define the couplings gy, g± as the following current-current perturbation of the su(2) level k = 1 Wess-Zumino- 
Witten (WZW) model: 

(ig ± (J+T + J-J + ) - 4 5 |, J 3 J 3 ) , (8) 

where J a (J ) are the left (right) -moving currents, normalized as in [|J. 
Let us first consider the well-known one-loop result. The beta functions are 

(») 

where the length scale is L = ae l , with a a microscopic distance. We will refer to I as the RG 'time'. The flows 
possess the RG invariant: 

Q = gf\- g±, (io) 

so that the RG trajectories are hyperbolas. We are interested in the region g\ > grjj where Q is negative. Let us 
parametrize Q as follows: 

Q = VQ^i~, (ii) 



where h > 0. Eliminating g±, one finds 



This is the same beta- function as in The latter is easily integrated: 

<?||(0 = -^tan(/i(Z-Z o )), (13) 

where is a constant. One sees from eq. I|13|l that gy flows to — oo then jumps to +oo and begins a new cycle. The 
periodicity of the RG flow is g\\(l + Ai_i oop ) = g\\(l) where Ai_i oop = w/h. 

Though the 1-loop RG already indicates a cyclicity, since the flows extend outside of the perturbative domain 
|<7|| I < 1 one cannot conclude the flows are indeed cyclic based on the 1-loop approximation. However it was shown 
in |j| that the cyclicity persists to higher orders in perturbation theory, as we now review. 

The all-orders beta function proposed in ,6l] is 

dg\\ _ -4gj(l + <?n) 2 dg± _ -ig±(g\\+gj) , 
dl (i-gir ' dl (l- 5 i)(l-.9||)- 1 j 
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A few remarks on the status of the conjectured beta-function eq. (|14(l are in order. A number of important checks 
were performed in £|. The most sensitive check was of the massless flows that arise in the imaginary sine-Gordon 
theory defined by g± — > ig± . The above beta-function correctly predicts the known non-perturbative relation between 
the anomalous dimensions in the UV and IR. An all orders beta-function was also proposed by Al. Zamolodchikov [jj, 
and his result was quoted in [lfj. Zamolodchikov's argument was global in nature and did not rely on summing up 
perturbation theory; the main input was the known properties of the massless flows, thus his beta-function appears 
to be the unique one up to a change of coordinates that captures the non-perturbative aspects of the flows. It can be 
shown that the beta-function in |l0(| is equivalent to l|14|) under a change of coordinates. Recently it has been argued 
that there are additional contributions that first arise at 4-loops ^1 ■ Though this would seem to spoil the properties 
of the massless flows, arguments were given in that the new contributions would not. This is a rather technical 
issue and consequently we will assume the above beta-function is correct. Clearly, the present work is a further check 
of the global properties of the conjectured all-orders beta function. 

The RG flows were analyzed in detail in 0]. It was shown there that though the beta function possesses poles 
at ±1, the flows approach the poles along non-singular trajectories and can be extended to all length scales. The 
analysis of the flows is simplified by recognizing that they possess a non-trivial RG invariant: 



Q 



9l 



(i + .9ii) 2 (i-si) 



(15) 



(We have rescaled Q by a factor of 16 in comparison to |4|.) 
One can use Q to eliminate g±, obtaining 



dg\\ 
dl 



= -4 



i-(i + .gn) 2 g) (^-(i + .gn) 2 Q 



(16) 



This is readily integrated 
1 



4VQ 



tanh 



g\\(l-Q)-Q 



— tanh 



Q{g\\ 



log 



/ .9 2 -(l + .9||) 2 Q S 
1-(1+.9||) 2 Q 



(17) 



where Iq is a constant. Note that there are no singularities at the location of the original poles of the beta function. 

As described in the RG flows possess a number of phases depending on the value of Q. For the flows with UV 
or IR fixed points the perturbing operators away from the fixed point are scaling fields, i.e. their scaling dimension is 
constant along the flow. Thus for the flows with fixed points, the RG invariant Q encodes these scaling dimensions. 
The flows with fixed points all correspond to Q > and are reviewed in more detail in the next section. 

The cyclic flows on the other hand correspond to Q < 0. Here, since there are no UV or IR fixed points, Q does 
not encode anomalous dimensions, but rather the only universal feature of the flow, namely the period of the cycles, 
which we now calculate. The coupling g» flows toward — oo, where it jumps to +00 and eventually returns to its initial 
value. One can see explicitly this jump from the solution as follows. Parameterizing Q as in eq. 1|11|) . for h large, and 
I on I w 00, the solution is 



■911 



h 



tan 



-Ml ~ h) 



(18) 



Thus when I — Iq — Tt/h the coupling jumps from —00 to +00. 
As before, define A as the period in I for one cycle, eq. QJ. 



Evidently X — l-oo — loo where l±ao are the RG times 



when gu = ±00. From the exact solution eq. i|17|) . and the fact that tanh 1 (ioo) = «7r/2, one finds the simple result 



A = 



2tt 
~h' 



(19) 



Note that this is twice the 1-loop result. 

The RG results presented so far also give an indication of the mass gaps of this theory. Indeed from eq. H18|l we see 
that when I — Iq = ir/h the coupling g\\ becomes —00, which implies a mass gap M ~ \ e ~' K ^ h associated to the length 
scale Be"/' 1 (we are assuming for simplicity Iq -C ir/h and h -C 1). This mass is usually associated with the existence 
of two particles, the spinons, in the application of this model to the XXZ spin chain in the antiferromagnetic regime 
(see section V). 
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III. BOSONIZATION 



The level-1 su(2) current algebra can be bosonized as follows: 



J = -=e 
V2 



±i\ r 2ip 



V2 



(20) 



where <p(z) is the z = t + ix-dependent part of a free massless scalar field cf> = ip(z) + tp{z). Viewing the g» coupling 
as a perturbation of the kinetic energy term and rescaling the held (j) one obtains the action 



(21) 



where to lowest order A oc g±. The held cf> is normalized such that when A = 0, (4>(z, ~z)<f)(0)) = — logzz and the scaling 
dimension of cos 60 is b 2 . (The coupling b is related to the conventional sine-Gordon coupling (3 as b 2 — /3 2 /4-7T so 
that the free fermion point corresponds to b = 1.) It was shown that the all-orders beta- function eq. H14fl is consistent 
with the known two- loop beta- function of the sine-Gordon theory in |12|. 

For the flows with fixed points, i.e. Q > 0, the coupling b can be related to Q by properly matching the slope of 
the beta function at the fixed points where g± = 0. The result is 3 



2b 2 



l + 2^Q 



b 2 ' 



(22) 



The region < \[Q < oo is the massive sine-Gordon phase with < b 2 < 2. If < b 2 < 1, i.e. y/Q > 1/2, the 
spectrum contains bound states of the solitons and antisolitons called breathers, which are absent if 1 < b 2 < 2, i.e. 
< \J~Q < 1/2. When y/Q < —1/2, b 2 < 0, and b is purely imaginary, this corresponds to the massive sinh-Gordon 
theory. Finally when — 1/2 < ^/Q < 0, the perturbation is irrelevant since b 2 > 2 and this corresponds to the massless 
regime with an IR fixed point. All these regions are plotted in figure 2. 



cyclic regime < 



sinh-Gordon 



-1/2 



ih/4 



1/2 



IR-massless 



no breathers 



breathers 



massive sine-Gordon 



FIG. 2: Description of the different regimes of the sine-Gordon model as a function of \[Q. 



Parameterizing Q as in eq. 
corresponds to 



(|1 1|) . and making an analytic extension of eq. I|22|l . we see that the cyclic regime 



b 2 = 



-itan _1 (/i/2) 



l + ih/2 (A + h 2 ) 1 / 2 



(23) 



Since b has both real and imaginary parts, this theory is neither the sine-Gordon nor the sinh-Gordon theories. In this 
bosonized description the hamiltonian does not appear to be hermitian however this is an artifact of the bosonized 
description since in the current algebra description the action is real. To lowest order, the origin of the imaginary 
part of b is the square-root in the rescaling of the held involved in going from eq. (|20|l to eq. (|21|l . Furthermore the 
spin-chain realization in section VI also has a hermitian hamiltonian. 
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IV. QUANTUM AFFINE SYMMETRY AND S-MATRICES 



In this section we will use some properties of the bosonized description to put forward a conjecture concerning the 
spectrum and S-matrices. The sine-Gordon theory is known to possess the U q (sl{2)) quantum affine symmetry which 
commutes with the hamiltonian and the S-matrix [13|. where 



2*i/b 2 



(24) 



The construction in [Tjj is valid regardless of value of b. For the cyclic regime, though b has real and imaginary parts, 
q is real: 



q = -e-^ 2 . 



(25) 



The form of b in eq. lll'-it is rather special in that it leads to the real value of q in eq. I|25|l : a more general complex b 
would imply a complex q. 

We point out that the sine-Gordon theory is also known to possess another, dual, quantum group structure Lh(sl(2)) 
which encodes the commutation relations of the monodromy matrix |l4lll5| . The dual quantum group is the algebraic 
structure underlying the Bethe ansatz in the quantum inverse scattering method. The latter is closely tied to the 
hamiltonian itself, since the trace of the monodromy matrix yields the integrals of motion. Generally, if q — e wa then 



the dual is q 



— p 1 /" 



Let us flip the sign of q, which in the field theory is nothing more than a Klein transformation 



modifying the statistics of the soliton fields. Then the dual-g is precisely the inverse of the finite RG transformation 
of one cycle: q~~ x = e 2 ""/' 1 = e x . Though this is an intriguing observation, we will make no further use of it in this 
paper. 

As for the usual sine-Gordon theory, we assume the theory contains solitons of topological c harg e ±1. In the XXZ 
spin chain, these solitons are the well-known spinons. (See the next section.) It was shown in [13j that the S-matrix 
is fixed up to an overall scalar factor by requiring it to commute with the quantum affine symmetry. We now describe 
the resulting soliton S-matrix in detail. Let (3 denote the usual rapidity that parameterizes the energy and momentum 
of the soliton: 



E = M s cosh 0, p = M s smb. (3. 



(26) 



Let A a ((3) denote the creation operators for the solitons, where a = ± denotes the topological charge. The S-matrix 
may be viewed as encoding the exchange relation of these operators: 



A*{fr)Ai{p2) = SlUh - ft) A d ((3 2 )AM. 
Requiring that the S-matrix commute with the quantum affine symmetry leads to the following structure: 



SI 



S(0) = 






V o 








-+ 








szz) 










where 



-nh/2 



_ p -ifih/2 



(27) 



(28) 



(29) 



([3 = Pi — 02-) The above S-matrix automatically satisfies the Yang-Baxter equation. 

The quantum affine symmetry does not constrain the overall scalar factor p, however there are additional constraints 
coming from crossing symmetry and unitarity. Crossing-symmetry requires: 



1 



p(C) = p(--) 



Unitarity corresponds to S j '((3)S(f3) = 1. If the S-matrix is real analytic, i.e. 

S^(0)=S(-(3) 

then unitary requires 



(30) 



(31) 



S(-(3)S((3) = 1 



(32) 
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The above equation leads to the following constraint on the overall scalar factor p: 

PiOpiC 1 ) = (1 _ cV) Ji_ c _ V) (33) 

A minimal solution to the above equations is easily found iteratively. (See e.g. H3-) One starts with a solution to 
eq. Ip!3^1 . pn = Hal (1 — a 2 C 2 \. This is made crossing symmetric as follows: p\ = po(C)i°o( — V^O- This spoils unitarity, 
which can be corrected by multiplying by pi(C _1 ). Then one starts over again correcting crossing symmetry. The 
process never terminates and leads to the infinite product: 

Kg ft (l-g 4+4 "C- 2 )(l-g 2+4n C 2 ) m] 
P[q, ^ ! 1 - q 2 C 2 W (1 - 9 4 + 4 ™C 2 )(l - g 2 + 4 "C" 2 ) 

This infinite product is convergent since |g| < 1. 

The above S-matrix is real analytic eq. (|31[) and thus unitary in the usual sense. Most remarkably, it possesses a 
periodicity that is a clear signature of a cyclic renormalization group. Up to some minus signs, the S-matrix satisfies: 

S(f3+ T ) = S((3) (35) 

In the limit of large energies, where E rj me^/2, the above equation implies eq. J3J. 

That the S-matrix presented in this section correctly describes the current-current perturbation in the cyclic regime 
is the main conjecture of this paper. The lagrangian defines a hermitian theory, and our S-matrix is unitary as it 
should be. It also has a periodicity that is correctly predicted by the RG analysis of section II. In the next section we 
provide further support of our conjecture based on the XXZ spin chain. 



V. REALIZATION IN THE XXZ HEISENBERG CHAIN 



In this section we give further support of the S-matrix proposed in the last section by taking a particular continuum 
limit of the XXZ spin chain. It is well-known that in the continuum limit the spin 1/2 Heisenberg chain is well 
described by the anisotropic current-current perturbation of section II. (See for instance 10, 30].) 

The spin chain has the hamiltonian 

N 

H = - J H KX+i + «+i + A ■ (36) 

i=l 

Let us first relate A to the couplings gn, g± in the critical (massless) regime, — 1 < A < 1. In the current-current 
RG description, here the flows terminate in the IR along the critical line < gn < 1 with g± = 0. In this regime, 
— 1/2 < \[Q < 0. The spin-chain may be mapped onto a theory of fermions via the Jordan- Wigner transformation. 
Bosonizing these Jordan- Wigner fermions, one obtains the sine-Gordon theory. (See for instance |lOLl30l l3i|.) In this 
bosonized description, the flows arrive at the fixed point via the irrelevant operator cos(6</)) with scaling dimension 
b 2 . The parameter b is known from the Bethe ansatz to be related to A as follows: A = cos27r/fe 2 . In terms of Q: 

A = -cos27tv/Q. (37) 

Let us turn now to the cyclic regime, where — oo < Q < 0. Let us assume that the formula 1)37(1 is valid in this 
regime also. Although ^/Q is imaginary, A is still real. Parameterizing Q as in eq. I|ll|) . ^/Q = ih/4, one finds 

A = -cosh^. (38) 

From eq. (|38|l one sees that the cyclic regime of the RG should correspond to the antiferromagnetic domain A < — 1 
in the spin chain. In terms of A, the period A of the RG flows is given in eq. (JHJ, namely A = n 2 / cosh~ 1 (— A). 

There is one additional check of the above series of mappings. The spin chain is known to directly possess the 
U q {sl{2)) symmetry on the lattice ^(|, where q is related to A as follows: 



(39) 
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Observe that the above equation is consistent with eqs. and l|2Bfl. 

The XXZ spin chain with A < — 1 was studied in using the underlying U q (sl{2)) symmetry on the lattice. One 
feature of this algebraic approach is the construction of so-called type II vertex operators Z a , a — ±1, which can be 
interpreted as creation operators for the fundamental spinons. These operators satisfy the exchange relation 

z a (Si)z b (&) = i$(&/G0 z d (( 2 )z c (Ci). (40) 

(Compare with eq. I|27(l ). The matrix R in the above equation can be interpreted as the S-matrix for the spinons on 
the lattice. An S-matrix for this regime of the spin chain was also obtained from the Bethe ansatz in |3^| . The explicit 
form of R is precisely as in eq. I|28fl where now q parametrizes A as in eq. (|39|l and is thus the same as in eq. Q25p. 
The expression for p on the lattice is the same as in eq. (|34|l ; in the algebraic construction of|l6j|. the S-matrix is just 
the universal R-matrix for the quantum affine algebra. 

On the lattice, the particles have energies expressed in terms of elliptic functions. Letting 

C = -*e M , (41) 

then in terms of a the energy and momentum are 

. . 2K . , irK' , (2Ka\ , . ( '2Ka\ it 
E(a) = sinh -^dn , p(a) = am - -. (42) 

7T K \ IT J V 7T J 2 

(See [33113] f° r definitions.) In the above formulas the nome of the elliptic functions is |g| = e~^ h l 2 = e _,rA l K which 
is between and 1. 

Let us take a particular continuum limit which amounts to letting the momentum and also h to be small. Let 
us first take the momentum to be small. From the definition of am(u), one finds that p(a = tt/2) = 0. Letting 
a — a — 7r/2, from the definition of am(u) one has 

K+ *KS_I~» # (43) 

which can be transformed to 

2Ka_f> d# 
K Jo \/k' 2 + k 2 sin 2 ip' 

Above, k, k! are the standard moduli of the elliptic functions. When p is small, sin^' can be approximated by ip' and 
one finds 

k' [2Kka\ . x 

p = — sinh . (45) 

k \ n J 

We now let \q\ approach one, i.e. h small but not zero. Using known expressions in terms of the dual nome 
q = e~ 27r /' 1 = e~ A , which goes to zero in the limit [34]], one has 

A-«£(l + 4g), K'tt^(l + 4q), fcwl-8§; k' w 4^. (46) 

Thus in this limit one finds 



where the rapidity /3 is identified as 



Turning to the energy, we need the identity 



-fc'sinh/3, (47) 



/3 = -f- (48) 



dn(u) = J 1 - k 2 sin 2 (am(u)). (49) 

Using the limits eq. I|46|l one finds 

E = Trfc'cosh/3. (50) 
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Thus the spinon velocity is it and the mass goes with k' = Ay/q = 4e _Tr / / \ The latter result agrees with the RG 
estimate of the mass gap Mo ~ ^e~ 7r /' i made in section II. Identifying the rapidity {3 from eq. H47|) one finds that the 
spectral parameter £ = e ~ 1 ^' 1 / 2 , which agrees with eq. (|29|) . 

To summarize this section, we have shown that a particular continuum limit of the XXZ spin chain exists which 
leads to the S-matrix conjectured in the last section. We also emphasize that in taking the particular continuum limit 
described above, we have implicitly provided an ultra-violet completion of the XXZ spin chain by defining a limit 
where the dispersion relation is relativistic. 



VI. ANALYTIC CONTINUATION OF THE SINE-GORDON S-MATRIX 



In this section we investigate another solution to the constraints in section IV on the S-matrix. Based on the 
bosonized description of section III, it is natural to consider the analytic extension of the sine-Gordon S-matrix to the 
complex value of the coupling b given in eq. We will refer to the resulting model as the cyclic sine-Gordon model. 

Since the usual sine-Gordon S-matrix also commutes with the quantum affine symmetry, the only difference between 
the cyclic sine-Gordon S-matrix and the one proposed in section IV is in the overall scalar factor p. Extending the 
formulas in |17| to the regime l|23[) . one finds the S-matrix is given by eq. (|28|l but with p given by psc, where: 

pOT» - -ir («) r<i - , )r (i - « + .) fl ^fgfe/ 2 . CD 



with 



where we have defined 



= T(2nn- z)T(l + 2nn- z) 

n[P> r((2n + l)«-*)r(l + (2n -!)«-«)' 1 ' 



ih ./3k (3h 

« = — , z = -% — = — . (53) 

Z 7T Z7T 

The standard S-matrix for the massive sine-Gordon model is given by the same formulas (|28I51I52|) with the following 
substitutions: 

«=^, z=-** : ?= _ e *rV7, C = e- S ^\ (54) 
7 7 

where 7 is related to b by the equation 



1 



in 2-b 2 2y/Q' 



VQ > 0. (55) 



Hence the relationship between the massive sine-Gordon S-matrix and the S-matrix considered in this section amounts 
simply to the analytic continuation 

h «-> (56) 
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The above S-matrix satisfies the constraints of section IV, including algebraic unitarity eq. (|32|l . however unlike 
the S-matrix in section IV it is not real analytic eq. (|3 1|) . and thus cannot be considered unitary in the usual sense. 
As discussed in |20j . in theories where parity is broken, real analyticity can be violated, and replaced by hermitian 
analyticity. For simplicity consider a diagonal scattering theory where the S-matrix for the scattering of particles 
A ai ((3i) and A a2 {(3 2 ) is given by S aia2 (fti - A-)- For a parity non-invariant theory, S aia2 (/3i - fo) ^ S a2ai ((3i - j3 2 )- 
Unitarity, S — 1, requires 

SZ aai (Pi-lh)S ai M-lh) = l (57) 
Imposing the requirement of hermitian analyticity 

S a2ai (Pi ~ fc) = S* aia2 ({3 2 - fc) (58) 
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the constrain of unitarity becomes 

Sa lQ2 (/3)S aiQ2 (-/3) = l (59) 

Note that the latter corresponds to algebraic unitarity eq. I|32|) which is satisfied by the cyclic sine-Gordon S-matrix. 
We will henceforth assume the lack of real analyticity can be understood in this way. 

The analytic extension eq. I|23[l is reminiscent of the staircase models however there are some important 

differences. The staircase model corresponds to b 2 — —2(1/2 + ih$ / 2ir) / (\ / 2 — ih^jl-K). In terms of \/Q, it corre- 
sponds to t/Q = —1/(1 + ih Q /ir). For small /i , this is a deformation of y/Q = —1, which is in the sinh-Gordon 
regime. The S-matrix is thus an analytic extension of the sinh-Gordon S-matrix, with a single resonance, whereas 
ours is a deformation of the sine-Gordon one. Furthermore, for the staircase model, the c-function, though it roams, 
monotomically decreases. For a model with a cyclic RG we expect a cyclic c-function (see the conclusion). 

A. Resonances 

The S-matrix for the cyclic sine-Gordon model has numerous poles which can be interpreted as resonances. Inter- 
estingly, other known models that violate parity real analyticity also have resonances [2£|. More generally, consider 
a pole at = \i — irj with /i, r\ > 0, in the S-matrix for the scattering of two particles of masses mi, m^- As discussed 
in 0, 122, H3| this corresponds to a resonance of mass M and inverse lifetime T where 

— — ^ = m\ + m 2 + 2miTO2 cosh(/x — irf). (60) 

Equivalently: 

r 2 

M 2 — = m\ + m\ + 2mim2 cosh/i cos 77, 

Mr = 2mirri2 sinh/isin?7. (61) 

The soliton S-matrix has an infinite number of poles at (3 n — 2ixnjh and n = iir — 2Tm/h with n a positive integer. 
These poles are absent in the S-matrix proposed in section IV. The poles at n yield a real mass and correspond to 
resonances with an infinitely long lifetime r = 1/r = 00, while the poles at n yield an imaginary mass and are not 
physical (see table 1). It is interesting to compare this situation to what happens for the usual massive regime of the 
sine-Gordon theory where < b 2 < 2. In this case, when 7 < 871", the poles at n lead to the breather bound states, 
while those at n are these same poles in the crossed channel. 



Cyclic sine- Gordon 


Massive sine- Gordon 


Resonances 

a 27m 
Pn — h 

m n = 2M S cosh ^ (n = 1, . . . , 00) 


Breathers in crossed channel 

a — in i 

Pn — 8 


Resonances in crossed channel 
n =iv-^ 


Breathers 

m„ = 2M s sin^ (n = l,...<^) 



Table 1.- Poles and masses of the cyclic and the massive sine-Gordon models. 

Since the poles n are right on the cut in the s-plane (s = (pi +P2) 2 ), it is desirable to incorporate a small imaginary 
part —irj taking them off the cut and giving the particles a large finite lifetime. This can be viewed as a kind of "ie" 
prescription as follows. Letting 

h^h + ie, (62) 

with e very small, the poles are now at 
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The prescription eq. I|62[) does not spoil crossing-symmetry, algebraic unitarity, nor the Yang-Baxter equation. As 
e — * 0, from eq. I|tjl|) the mass m n and width T n of the n-th resonance is 

m n 2M S cosh — , r„ w e ~^r M s sinh — . (64) 

Each resonance of mass m„, n = 1, 2, ..oo, actually corresponds to four particles transforming in the direct sum of 
the q-deformed spin 1 and singlet representations. One can easily see this by noting that the entries S++,Sl[_ and 
their charge conjugates are non-zero at the poles /3 n , while Stz and SZ+ are zero. If the fundamental solitons have 
topological charge ±1, then each resonance is a triplet with topological charges (—2, 0, 2) or a singlet with topological 
charge 0. This is in contrast to the usual breather bound state poles at f3 n which occur in the scattering channel of a 
soliton with an anti-soliton and are thus singlets. (Sff and SZZ are zero at the pole (3 n ). 

In the language of the fusion procedure |2^|, which is equivalent to the bootstrap, the poles j3 n project onto the 
singlet whereas (5 n projects onto the (g-deformed) spin-1 and spin-0 representations. This can be seen explicitly by 
going to the so-called homogeneous gradation. The soliton S-matrix can be expressed as 

S(C,q)=-^-.<J2iR<J u 1 P, (65) 

H /2 H /2 H /2 H /2 

where on = Ci ® C2 > a 2i = (2 ® Ci 1 with H = diag(l, —1) the topological charge and P the permutation 
operator: 



P= • (66) 




The matrix R can then be expressed in terms of the projectors Pq and Pi onto the q-deformed spin and spin 1 
representations : 

R(0 = ((q - CV 1 )^ + (C'q - (q-^Po, (67) 

where 

/0 s 
1 -q 



l + q 2 



-q q 2 j ' (68) 

Vo oo, 



and Pi = 1 — Po. The usual sine-Gordon breather poles at n correspond to ( — (— )" +1 (7 _1 and one sees that at 
this C, R oc Pq, showing that the breathers are indeed singlets. A pole at £ = q projects onto the irreducible spin 1 
representation. However our resonance poles at /3 n correspond to £ = (— ) n where R oc Pq + P%. The resonances, as 
the sine-Gordon breathers 0]> also carry a C-parity given by (—1)" for the triplets and (— 1)™ +1 for the singlets. 

The above spectrum of masses has scaling properties that are consistent with a cyclic RG. As argued in section III 
the spinon mass M s must be proportional to l/L, where L is the system size. In the limit n 3> h/n or h <S 1, eq. 164|) 
becomes, 

m„(£)~)e M/ '', n^h/n, (69) 

which exhibits the Russian doll scaling property 

m n {e- x L) m m n+2 (L), (70) 

similar to what was found in Thus the way the spectrum reproduces itself after one RG cycle is that m n at a 
given length e~ A L plays the same role as m„+2 at the longer length L (see figure 3). Notice that the jump by two 
in n is consistent with the C— parity of the resonances. Eq. (j70|l also shows that after one RG cycle two new low 
energy masses appear in the spectrum. This gain of the lowest energy states in the spectrum is what allows for the 
reshuffling of resonances after one cycle and seems to be an essential ingredient of any cyclic RG. 

In a finite size system we expect to have a finite number of resonances, n rcs , which can be estimated as follows. 
If in every RG cycle we get two more resonances, then after p cycles with L ~ ae pX , we expect to have n rcs = 2p 
resonances, with 

h, L 

"-res Wg-. (71) 

7r a 

This kind of equation was checked numerically for the Russian Doll BCS model[f|. 
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FIG. 3: Log-log plot of the resonance masses m n (L) versus 1/L, where L is the length scale of the system. Lo denotes the 
maximal size of the system. The Russian doll scaling of eq. (I7UH is described by the dotted lines. 



B. Closing the bootstrap: stringy spectrum 

In the limit that the regulator e goes to zero, the resonances are stable particles and can exist as asymptotic states. 
Thus, one should attempt to close the bootstrap including the resonances in asymptotic states as proposed in |23|. 
Since the bootstrap is equivalent to the fusion construction in the theory of quantum affine algebras, it is clear for 
example that the resulting S-matrix for the scattering of the soliton with the triplets of resonances should be solution 
of the Yang-Baxter equation corresponding to the universal 7?.-matrix of U q (sl{2)) evaluated in the representations 
V1/2 8> Vj. where Vj is the 2j + 1 dimensional g-deformed spin-j representation. 

Continuing to close the bootstrap leads to resonances in the q-deformed spin j representation Vj with j = 
1, 3/2, 2, ...oo. To simplify matters we will only consider the scattering of the particle in Vj of maximal topologi- 
cal charge +2j. This avoids the complexities of non-diagonal scattering and will be sufficient to obtain a mass formula 
as a function of j. 

Let 5"j™2 j W) denote the diagonal scattering of the 7i-th spin j particle of topological charge 2j with the soliton 
of topological charge +1. For j = 1, n refers to the quantum number in eq. I|64|) . and when j = 1/2, there is no 
n-dependence and Sy 2 t , 2 = If the particle of charge 2(j + 1/2) appears as a resonance in the scattering of spin 

1/2 with spin j particles, then the S-matrix has a simple pole: 

S i/l,j(P) ~ ~ 1+172' ( 72 ) 

P Pl/2,i 

and the masses are related by the formula 

m 2 c = m a + ml + 2m a mb cosh fi^ b , (73) 
with a,b,c = 1/2, 1, 3/2, .... The above equation implies the following relation between the fusion angles: 

Vab + Mbc + Mac - 2™. (74) 

In the equations 173I74|) all the masses are implicitly at the same fixed n. 

Given S^l j one can construct the S-matrix S^ 2 j +1 / 2 fr° m the bootstrap procedure: 

S l/2 .j + 1/2^1 3 ) = S l/l,j + PjJ + l/2) 5 'l/2,l/2 [0 - / I l/2 J + l/ 2 ) ' ( 75 ) 

where ~p — iir — \i. 

Let us first consider the scattering of the spin 1/2 solitons with the spin 1 resonances discussed above. We already 
found that ^1/2 1/2 = 27rn//i and eq. I|74|l implies 7*1/2 1 = Kn/h. The bootstrap equation then reads 

= Si/2,1/2 (p + ^) S 1/2A/2 (0 - . (76) 
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The factor Sx/2,\/2{P) om y has real poles at z = f3h/2n = p and zeros at z = —p where p = 1,2,3, ... Remarkably 
this leads to the fact that the only real pole of Sy^ is at f3 — nn/h = y^j^ v The other poles of Si/2,1/2 either 
lead to double poles or are canceled by the zeros in the other Si/2,1/2 factor. 

Let m n (j) denote the mass of the n-th resonance with spin j. The above value of M1/2 1 l e ads to 

ml(j = 3/2) = 4Af 2 ( 9 - + 2 sinh 2 ^ . (77) 



Moving on to the spin 2 resonances, one finds that the fusion angles are somewhat more complicated. Since the 
masses m(j) are already determined for j = 1/2, 1, 3/2, the relevant fusion angles follow from eq. (|74(l and one obtains 

3/ 2 _ mr 

^1/2,1 ~~ L ' 



««P^iA3/2 = V2T^' H = e 2 ™/\ (78) 



_ 1/2 _ / (l + ff)(2 + ff-i) 
exp Ml3/2 - W (1 + ff _i )(2 + ff) - 



The S-matrix is 



Similarly to the j — 3/2 case, the above S-matrix has a single simple pole at (3 = yf^ 2 3/2 where 



Mi/2,3/2 ~ lh%2- ( 80 ) 
This leads to the mass 

m 2 (j = 2) = 4M 2 (4 + 3 sinh 2 ^) . (81) 

Having understood the above cases we can readily extend the results to arbitrary j. Examining the poles in the 
S-matrix <Si/2,j one finds a single, real, simple pole at 

3+1/2 _ j _ 1/2 / R9 s 

Ml/2j -Ml/2,j-l/2 ^3-1/2,0- 

This leads to the mass formula 

m 2 (j) = 4M 2 (j 2 + (2j - 1) sinh 2 ^) (83) 
= M 5 2 (2j-l + ff)(2i-l + J H'- 1 ). 
From this mass formula one can determine all the fusion angles: 



j+1/2 2J-1 + H 

exp ^/ 2 , - V aj-i + g- 1 ' 



-7 2j + H 

cx P/4/ 2j+ i/2 = \/ 2j +g -i ' (84) 



_i/2 / (2j-l + g)(2j + g-i) 

ex PM JJ+ i/2 - y(2j-l + ff-i)(2j + H)' 

The above mass spectrum is suggestive of a string theory. We remark that the higher spin particles that make up the 
Regge trajectory of the bosonic string also appear as resonances on the physical cut in the Veneziano amplitude |25j | . 
Consider the "leading Regge trajectory", i.e. the lowest mass at spin j. When h is very small, the formula 183|) 
becomes 

m 2 (j)«^(j-a(0)), (85) 
a' 
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where the Regge slope and intercept are 

= 2M s 2 e 2 ^ ft , a(0) = 1/2. (86) 

cr 

(The leading Regge trajectory corresponds to n = 1.) The usual field theory limit in string theory is a' — > which 
corresponds to h — > 0, where all the resonances become infinitely massive. The fact that our theory has two coupling 
constants, namely M s and h, also suggests that the underlying string has another parameter such as a compactification 
radius. If there is indeed an underlying string description, much remains to be understood regarding the nature of 
this string theory. 



C. Jumping the c m = 1 barrier? 

We close this section with some speculative remarks regarding the possible nature of the string theory, if there is 
indeed such a description. Consider a non-critical bosonic string theory in D spacetime dimensions, where c m = D — 1 
is the Virasoro central charge of the matter content with c m < 26. Because c m ^ 26, there are some Liouville degrees 
of freedom with central charge cl, where c m + = 26. The Liouville theory has the action 

liouville = / ^ ( J W) 2 + \e* b A . (87) 



4tt V2 

The central charge of the Liouville theory is related to b as follows: 



ib V / 2 X 



Ci = l + 6^= + ^J . (88) 

The full range 1 < c m < 25 beyond the "c,„ = 1" barrier, corresponds to 

b 2 = 2e la , with cosa= Cm ~ 13 (89) 

12 v ' 

Note that in the range we are considering, —1 < cos a < 1. Like the sine-Gordon theory, the Liouville theory has a 
finite U q (sl(2)) quantum group symmetry with q the same as in eq. (|24() 13, 26] . 

To relate our S-matrix theory to the string theory we would need to relate the Polyakov path integral over Riemann 
surfaces (summed over genus) for the matter coupled to 2d gravity, and it isn't at all clear how to even begin to do 
this. We point out however that this has been done for the c m = 1 case leading to an exact 2D S-matrix |27J. (For 
a review, see|2^,|2^.) Here we merely examine the following. The e lh ^ interaction of the Liouville theory is one part 
of the cos b(j) interaction of the sine-Gordon theory, and we expect that the Liouville theory for b given in eq. H89|) 
shares some common features with the sine- Gordon theory at the same b. 

For a = —h/2, and h small, eq. (|89|l matches eq. (|23|) and thus corresponds to our cyclic sine-Gordon model. One 
sees that this is just below the c m = 25 barrier: 

Cm « 25 - ^h 2 (90) 

Interestingly, when a m ir — h/2 with h small, c m is just above 1, c m ~ 1 + 3h 2 /2, and this corresponds to the 
deformation of the sinh-Gordon theory which is the staircase model |18j- It would be very interesting if the idea of a 
Russian doll RG plays a role in resolving the known long-standing difficulties in going beyond c m = 1. 



D. Comparison between the cyclic sine-Gordon model and Russian doll BCS model 

It is interesting to make a comparison between the cyclic sine-Gordon model we have studied in the previous sections 
and the Russian doll BCS model introduced in reference 0, which also possesses a cyclic RG. The model considered 
in the latter reference is a simple modification of the reduced BCS model used in the study of the superconducting 
properties of ultrasmall metallic grains. It describes the pairing interactions between N/2 pairs of electrons occupying 
iV energy levels £j, which are separated by an energy distance 25, i.e. Sj+\ — Sj = 25. The model is characterized by 
a scattering potential Vjj' equal to g + i h for Sj > By and g — ih for Ej < eji in units of the energy spacing 5. 
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The case h = is equivalent to the usual BCS model which has a unique condensate with a superconducting gap 
given by A ~ 2N5e~ 1 / g for g sufficiently small. When h is non zero the BCS gap equation has an infinite number of 
solutions corresponding to condensates with gaps A n (N) ~ N A5e~ n7 '/ h (n — 0,1, ... , oo), where A depends on g and 
h. This model can also be studied using a renormalization group which reduces the number of energy levels N, namely 
N(s) — e~ s N, by integrating out the high energy modes. The RG leaves invariant the coupling constant h, while g(s) 
runs with an equation similar to the one loop result l|12|) and has a cyclic behavior with a period Xbcs — ^ /h, which 
in turn is related to the cyclicity of the gaps A n (e~ XBCS N) w A„ + i(iV). Notice that in this case, as compared to 
eq. (|70|) . the jump in n is one, due to the fact that there is no discrete symmetry associated to n. In a finite system 
the number of the condensates n c is given by n c ~ — logiV, in close resemblance to eq. H71jl. All these results together 
with the corresponding ones for the cyclic sine-Gordon model are summarized in table 2. 





Cyclic sine- Gordon 


Russian doll BCS 


RG-time 


L = e l a 


N(s) = e- s N 


RG-period 


ACSG = — 


Xbcs = f 




Resonances 


Condensates 


Energy scales 


m n (L) ~ y™/ h 


A„(iV) ~ NAe~ n7T / h 


Russian doll scaling 


m n {e' XcsG L) « m n+2 {L) 


A„(e~ ABCS N) « A n+ i(JV) 


Finite systems 


"res ~ £ log(L/a) 


n c ~ ^ log TV 



Table 2.- Comparison between the cyclic sine-Gordon model and the Russian doll BCS models. 

VII. CONCLUSIONS 

In this paper we have analyzed in detail the cyclic regime of the Kosterlitz-Thouless flows of a current-current 
perturbation of the su(2) level k — 1 WZW model [j]. We have computed the RG time A of a complete cycle in terms 
of the RG invariant Q = —h 2 /16 of the non perturbative RG equations Q, obtaining A = 2w/h. Using standard 
bosonization techniques we have mapped the perturbed WZW model into the sine-Gordon model parametrized by a 
complex coupling b 2 = 2/(l+ih/2). The latter model possesses a quantum affine algebra U q (si (2)) with a real quantum 
deformation parameter q = e~ vh l 2 , together with a dual quantum group structure U~(sl(2)) with q^ 1 — e 2 *"/' 1 = e A . 
Assuming the existence of solitons with topological charge ±1 with a factorized scattering S-matrix commuting with 
the U q (sl(2)) symmetry, we have examined two different S-matrix solutions differing by overall scalar factors. 

The first S-matrix we considered in section IV is real analytic and thus strictly unitary in the usual sense. It can 
also be obtained from a certain limit of the XXZ spn chain. The S-matrix possesses a cyclicity as a function of energy 
eq. which is a clear signature of a cyclic RG. We conjectured that it describes the cyclic regime of the KT flows. 

The other S-matrix we considered in section VI is an analytic extension of the usual massive sine-Gordon S-matrix. 
In this case the cyclic S-matrix has poles leading to an infinite number of resonances with the Russian doll scaling 
behavior also consistent with a cyclic RG flow. However since this S-matrix is not unitary in the usual sense because 
of the lack of real analyticity, it cannot describe the hermitian theory defined by the current-current perturbation. 
Thus it's physical meaning remains unclear. If there is a physically sensible underlying model, it probably corresponds 
to a parity violating theory as discussed in[2(j. Perhaps such a theory can be obtained as a certain limit of the XYZ 
chain. 

Remarkably we found that closure of the cyclic sine-Gordon S-matrix bootstrap leads to resonances of arbitrarily 
higher spin with a string-like mass formula. As discussed in section VI, it would be very interesting if the idea of a 
cyclic RG could shed new light on the problem of the c = 1 barrier. 

A few remarks concerning the so-called c-theorem are warranted. The c-theorem is a statement regarding the 
irreversibility of the RG flow: as the scale is increased massive states decouple and c decreases monotonically 35]. 
The c-function may be viewed as a function of the scale-dependent couplings, c(L) = c(g(L)). Thus, a cyclic RG 
would appear to give rise to a periodic c-function, violating the c-theorem. A thermodynamic Bethe ansatz analysis 
|3fij | of the c-function based on the S-matrix proposed in section IV indeed reveals a periodic c-function in the deep 
ultraviolet [37j. How can this violation of the fundamental c-theorem occur in a unitary theory as commonplace as 
the KT flows? One answer is that Zamolodchikov's proof of the c-theorem assumes that one can define the theory as 
a perturbation about a UV fixed point, however a cyclic RG by definition has no such fixed point. 

This work further supports the idea that cyclic RG behavior may be more commonplace than originally anticipated, 
and may possibly represent a new paradigm for Physics. It is tempting, though probably foolhardy to even speculate 
that the generations of fermions in the standard model are just the first few in an infinite sequence of Russian dolls. 
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